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SYMBOLS 


Coefficient defined by Equations (3) and (4) 
Random phase angle, 0 < B.S 27 


Coefficients in the equation x = Cl + C2 t +C3 rc used to 
search for a peak or valley value 


Harmonic coefficient, see Figure 1 

HARM subroutine control, IFSET 

= 1 sets up sine and inverse tables and calculates Fourier 
transform 


Integer representing time in HARM subroutine 


Mean value used in mapping Gaussian to simulated real time 
history 


N = oo see N 
Number of points to be computed in the random time history 


The number of points in the linear interpolation of subrou- 
tine TLU 


Factor by which the PSD frequency axis is extended to increase 
the sensitivity of the output 


Number of input data points on the PSD vs frequency graph 
Number of divisions of the PSD graph used in analysis 
Theoretical number of mean level crossings per unit time 
Theoretical number of peaks per unit time 


Theoretical number of Gaussian peaks exceeding the value x 
in the total time represented by the Gaussian time history 


Theoretical number of simulated real peaks exceeding the value 
y in the total time represented by the Gaussian time history 


Theoretical number of Gaussian peaks exceeding the value x 
per unit time 


Theoretical number of simulated real peaks exceeding the value 
y per unit time 


Probability of not exceeding x, normal probability distribution 


In order for a peak and valley to be counted and included in 
the output list, the maximum rise or fall must exceed the input 
value of PMAX, and the minimum rise and fall must exceed the 
input value of PMIN 


Vi 


PSD 


Power Spectral Density 


Transformation exponent used in mapping Gaussian to simulated 
real time history 


Difference in value of two consecutive Gaussian time history 
points, i, itl 


Difference in value of next two consecutive points, i+ 1, 
i+ 2, see SLI 


Power spectral density 


Time 


An independent value, also a Gaussian history value 


Independent variable to ue used in interpolation in a two- 
dimensional table 


A dependent value, also a simulated real time history value 


Dependent variable determined by interpolation of a two- 
dimensional table. 


Standard deviation 
Frequency 


Upper limit of frequency of input power spectral density 


Vii 


1. SUMMARY 


Load factor time histories developed using random noise simulation have 
been shown to correlate well with measured flight load factor histories. This 
makes random noise theory attractive for generation of fatigue test spectra. 

A common fault of other methods of load factor history simulation is arbitrary, 
unrealistic, coupling of peaks and valleys. The random noise theory approach 
requires that the random load factor time history be generated possessing a 
specified power spectral density (PSD) shape, and mean and RMS level. An 
advantage of the random noise theory approach is that both the exceedance con- 
tent and the frequency content.of the process can be preserved. The preserva- 
tion of the frequency content assures that proper coupling of peaks and valleys 
is attained. 

Volume II of this report, Reference 1, describes a computer program that 
uses as input lists of filtered peaks and valleys generated by this progran, 
and performs operations to combine and create variations of the input lists. 
Both of these computer programs were used in a study of load sequence effects 
on crack growth, summarized in Reference 2. The input data for the spectra 
generated for that study are listed in Appendix C, 

Gaussian Time History Generation 

The theory used in formulating the program is described in Reference 2. 

The basic equation used in forming the random stress histwry is, 
N-1 


x, = v2 2 cos (27u,t + BY) (2) 


where the term De is defined in Figure 1. Making the substitution De. = ¥S Aw 
and assuming Aw is constant throughout the summation described by Equation l, 


that equation can be re-written, 


N-1 
x, = ¥2Au Ls, cos (27w, + B.) (2) 
K=0 


Rather than solving this equation directly, considerable efficiency is attained 
by use of the Fast Fourier Transform (FFI) outlined in Reference 3. The 

FFT computes N points in the random time history for N input values of Avs 
where A. is computed, 


Re A, = 72Aw vs, cos B, (3) 


(J, = Upper Limit 
at Frequency 
of Input PSD 


S - Power Spectral Density - Stress7/cps 


Ww - Frequency, cps 


FIGURE 1. POWER SPECTRAL DENSITY USED IN 
RANDOM STRESS HISTORY GENERATION 


Im A, = Y2dw VS. sin By (4) 
Summarizing the generation of Gaussian random time histories: 
o "N' values of PSD (S.) are determined from the PSD curve corresponding 
£G- Kom: 0, Ly wag NEL 
o "N" random phase angles (B.) between O and 27 are generated. 
o An array of N complex coefficients, Ayes are calculated as shown in 
Equations (3) and (4). 
o The Fast Fourier Transform is used to generate a time history of N 
values of load factor x» as defined by Equation (2). 
Mapping of Gaussian Process to Real 
The exceedance curve for a Gaussian random time history is symmetrical 
with zero mean, as depicted in Figure 2. Actual exceedance curves for fighter 
aircraft are asymmetrical with a non-zero mean, also depicted in Figure 2, 
Gaussian random histories are adjusted so that exceedance curves for the Gaussian 
process and the actual exceedance curves are matched. This is accomplished 


using the transformations: 


y = (x + my ® for positive x + M (5) 


y= -(-x - my ® for negative x + M (6) 


Numbered points are used toe transform the 
Gaussian data to the sctual process 


N,, Exceedances - 1000 Hrs 


0 
X- % Design Limit Stress 


Exceedance Curve for Gaussian Random | 
Time History is Symmetrical with a Zero Mean 


Y = —(—x—mM)R 
for Y<O 


Ny. Exceedances - 1000 Hrs 


0 
Y - % Design Limit Stress 


Actual Exceedance Curves are 
Asymmetrical with a Non-Zero Mean 


FIGURE 2. COMPARISON OF EXCEEDANCE CURVES 


and conversely, 


x = y 


x 


eee M for positive y (7) 


~ (ay) B/R_ M for negative y (8) 


The exceedance curve for the Gaussian process is: 


nx” /267 


Ne = N,, P(x/ac) + Nf - P(x/bo) Je (9) 


This equation also describes the exceedance curve for the actual process when 
the substitutions indicated by Equations (7) and (8) are performed. The 
values M, R, N_ and o are unknowns, and are determined through a trial and error 
process by the user. Using the actual exceedance curves, two points on the 
valley curve and two points on the peak curve are selected. This results in 
four sets of values of N and y. By using Equation (9) in combination with 
Equations (7) and (8), four non-linear simultaneous equations can be written, 
relating N and y with M, R, ne and o as the unknowns. 

The ratio NA/N can be determined for any given PSD shape using Equations 


(10) and (11), 
f° 2 
w §(w) dw 


ppc gy hea Bae bees (10) 
. he S Cw) dw 
fe) 
Ou 
2 f i eCard 
N = ae Gea, TE (11) 
P f uw S(w)dw 
O 


The values of M and R are required for program input. The area of the PSD vs 


frequency plot is 9”, 5 ry ee 


u 
oo = fsa (12) 


The value of o can be employed by the user to scale the input values of PSD 


to obtain the proper area. 


search for Peaks and Valleys 


The time history generated by the technique outlined in the previous sections 


results in a series of stress values at constant time intervals, as de- 
picted in the following sketch. The presence of a peak or valley is detected 
when there is a difference in sign of slope between two consecutive 


sets of points. The value of the peak or valley is determined by assuming a 


quadratic relation of the forn, 


x= Cl+ C2 t+ C3 be? (13) 


and solving for the maximum value of x. 


Slope = SL, Slope = SL4 
Stress .* 
% 


1 


AT = Constant 


Rise and Fall Counting 


Counting accelerometers and VGH recorders are the primary basis of mea- 
suring the service experience of fighter aircraft, summarized in exceedance 
tables. This data, obtained from flight measured accelerations, has a rise 
and fall counting criterion applied to eliminate small load factor excursions. 
For example, VGH recorders require from .4g to 1.2g rise or fall, depending 
on the g level, before the peak or valley is counted. In a similar manner, 
the simulated real time stress time history is filtered, using input rise 
and fall criteria, in order to eliminate excursions that the user does not 


wish to consider. 


2. PROGRAM OUTLINE 


This program, based on random noise theory, generates random time histories 
for baseline spectra. The program is based on the theory described in 
Reference 2, and summarized in Section 1. Input to the program includes control 
numbers, the number of time points to be computed, the PSD as a function of 
frequency, an initial number to be used in the random number generator, the 
rise and fall criteria to be used for cycle counting, and Gaussian to simulated 
real transformation constants. Output includes a sequential list of peaks 
and valleys. Input and output are described in Sections 3 and 4, respectively. 

The procedure used in developing the stress random time history is out- 
lined in Figure 3. The first step is the reading of the input by the main 
program, and printing of the input. This is followed by computation of the A(K) 
coefficients described by Equations (3) and (4), where S(K) is the PSD, and 
B(K) is a random phase angle, computed using the techniques described in 
Reference 1. With the A(K) coefficients determined, subroutine HARM is called, 
which transforms the A(K) values into the Gaussian random time history. The 
main program stores the real part of the Gaussian time history in the last half 
of the A(K) matrix, searches for peaks and valleys and stores the peaks and 
valleys in the first half of the A(K) matrix and equations (9), (10), amd (11). 
The main program then prints the Gaussian occurrence and exceedance tables. 

Next, the main program tabulates the actual Gaussian peak and valley 
occurrence and exceedance arrays as determined from the list of peaks and 
valleys. Subroutine EXCED is then called, which computes the theoretical ex- 
ceedance array, based on the input PSD values. The main program then prints 
the Gaussian occurrence and exceedance tables. 

The next process is the transformation of the Gaussian process peak and 
valley stress history into a simulated real peak and valley stress history using 
the technique outlined in Section 1. The main program then tabulates the actual 
transformed peak and valley occurrence and exceedance arrays as determined from 
the list of transformed peaks and valleys. Subroutine EXCEED is again called, 
to compute the theoretical exceedance table, based on the input PSD values 
and transformation constants. The main program then prints the transformed 
peak and valley occurrence and exceedance tables. 

Subroutine FILTER is called to apply the input rise and fall criteria to 


the transformed, i.e., simulated real peak and valley stress history. The main 
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program then tabulates the actual occurrence and exceedance arrays based on the 
filtered transformed history. Subroutine EXCEED is again called, to compute 
the theoretical exceedance table, based on the input PSD values and transfor- 
mation constants. The main program then prints the filtered transformed peak 
and valley occurrence and exceedance tables. 

Finally, the filtered simulated loads are recorded on magnetic tape or 
disk for use by the program described in Volume II and then printed. 

The program consists of a main program and five subroutines, outlined 
in Figure 4. The function of each subroutine is described in the following 


paragraphs. 


Main Program 


Read Input Data 
Print Input Data 
Compute A(K) Terms 
Transform Gaussian History 
to Simulated Real 

Search for Peaks and Valleys 
Print Output Data 


Subroutine Exced Subroutine Harm 


® Compute Theoretical 
Exceedance Tables 
Subroutine NDTR 
@ Compute Normal! 
Probability Distribution 
Subroutine TLU 


@® Perform Linear 
Interpolation 


FIGURE 4, SUBROUTINES USED IN DEVELOPING STRESS 
RANDOM TIME HISTORY 


@ Transform A(K) Terms _ 


to Gaussian Time History 


Subroutine Filter 


@ Filter Peaks and 
Valleys 


Main Program 
The main program reads and prints the input data. The A(K) coefficients 


are computed, and subroutine HARM is called, which transforms the A(K) values 
into the Gaussian random time history. The time history is searched to deter- 
mine peaks and valleys. The main program transforms the Gaussian peak and 
valley values to simulated real values. The actual occurrence and exceedance 
arrays for the Gaussian values, transformed values, and filtered transformed 
values are tabulated by the main program. All output is printed or written by 
the main program. 
Subroutine EXCED 

This subroutine computes theoretical exceedance tables for the Gaussian 
and transformed, i.e., simulated real stress time histories. The solution for 
the Gaussian time history is obtained from the equation for exceedance per 
unit time, Equation (9), where Ny and N are computed from Equations (10) and 


(11). The Gaussian exceedance is computed 


Mota lice ~ . ss ve (14) 


where t is the time period represented by the input values of PSD maximum 
frequency wy and the number of computed time points. The theoretical exceed- 
ances for the transformed stress history is also obtained from Equation (9), 
with transformation indicated by Equations (7) and (8), 


Nrotal,y = tx N (15) 


Subroutine TLU 

This subroutine provides a linear interpolation table look-up procedure. 
Input to the routine is a list of independent and dependent values, x and y, 
the number of points in this list, NOPTS, and the independent value for which 
interpolation is desired, XVAL. Output of the routine is interpolated value 
of the dependent value, YVAL. When the range of the table is exceeded, 
linear extrapolation is used. 
Subroutine FILTER 

This subroutine applies the rise and fall criteria to the transformed 


peak and valley list. 


Filtering is accomplished by using the input values P MAX and P MIN on 
adjacent half cycles throughout the time history. A peak and valley that 
meets either the P MAX or PIN rise and fall criteria is counted only if the 
succeeding half cycle exceeds the remaining rise and fall criterion. If a half 
cycle fails to exceed the proper rise or fall requirement, the non-qualifying 
peak or valley is removed from the time history. The two adjacent peaks or 
valleys are now compared, resulting in the retention of either the largest peak 
or the smallest valley. The process continues by checking to see if the next 
half cycle exceeds the remaining rise and fall criterion. For example, in Figure 
3, ranges 4~5, 6-7, 7-8, 9-10, 10-11, 11-12, and 12-13 do not meet the rise 
and fall criterion. Sequential applications of the filtering checks results 
in the load sequence shown as "Fourth Pass". 
Subroutine HARM 

For one dimension, as used in this program, the HARM subroutine, described 


fully in Reference 3, computes the complex Fourier series, 


N-1 
KJ) = SA(K) oem JK/N (16) 
K=0o 


where J = 0, 1 ~---- N-1 


The theory of the transformation is described in Reference 3. Input to the 
subroutine consists of N, the real and imaginary parts of A(K), and the value 
of IFSET. Although the routine is capable of transforming three-dimensional arrays, in 
this program only a single dimensional array is required. Therefore, 
M(1) = Ml where N = gil is the number of time points to be computed, and 
M(2) = M(3) = 0. The real and imaginary parts of A(K) are computed by the 
main program, as previously described. IFSET is an option parameter which is 
set to 1 in order to have the routine calculate the Fourier transform. 

Output of the routine is the history X(J) stored in the A(K) matrix. 
The real part of A(K) is stored in even numbered locations and the imaginary 
part of A(K) is stored in odd numbered locations, sequentially. Subroutine 
HARM is a portion of the IBM Scientific Subroutine Package. 
Subroutine NDTR 

This subroutine is called by subroutine EXCED, in order to compute the 


theoretical exceedance table. The subroutine computes the normal probability 


distribution 
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Unfiltered 


First Pass 


Second Pass 


Third Pass 


Fourth Pass 
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FIGURE 5. APPLICATION OF FILTERING 


RULES TO STRESS HISTORY 


— 2 
y =. P(x) = = exp (-u' /2)du (17) 
as 


The following approximation is used: 


P (x) = 1 - £(x) > a,w's x > 0 (18) 
i=l 
where 

W = 1/(1 + px) 

f(x) = exp (-x*/2)/v21 

p = 0.2316419 

a, = 0.3193815 

a, = -0.3565638 

a, = 1.781478 

a -1.821256 

as = 1.330274 


The theory is described in Reference 4, Subroutine NDTR is a portion of the 
IBM Scientific Subroutine Package. 
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3. INPUT DATA REQUIREMENTS 


Input to the program consists of the number of time points to be computed, 
the power spectral density (PSD) as a function of frequency, the rise and 
fall criteria to be used with cycle counting, and an initial number to be used 
in the random number generator. The input data format is described in the 
computer program listing, via comment cards, Appendix A. Example input is 
presented in Appendix B. 

The maximum number of time points that may be computed in a single case 
is ra = 32,678. The number of time points to be computed is controlled by 
an input control number, and may be any power of 2, up to the maximum value. 

The power spectral density is usually input as (% Design limit weress)* 
vs frequency, cps. Program input includes the number of points to be input, 
and the frequency and PSD value for each of the input points. Experience has 
shown that many peaks and valleys will be missed if the upper limit of the 
PSD plot is not artificially extended with zero values for PSD, as depicted in 
Figure 6. The time represented by the input values of PSD upper frequency and 


the number of time points to be computed is, 


By increasing w artificially, the time t is reduced for a given number of time 
points, N, effectively increasing the number of points per unit time and hence 
per peak. Trial runs have shown that the use of NR' = 7 = 4, as defined in 
Figure 6, will usually result in adequate sensitivity and peaks and valleys will 
not be missed. 

The HARM subroutine requires as input a value of A(K) for each time point 
computed. Each value of A(K) requires a value of PSD to be interpolated. 
Because the number of time points to be computed, typically 32,768, is usually 
much larger than required to adequately define the shape of the PSD plot, 
constant values of PSD are used in ranges, as shown in the example of Figure 6. 
In that example, the number of time points is 128 and the frequency scale has 
been multiplied by four to increase the sensitivity of the output. Therefore, 
128/4 = 32 values of PSD are required in the computation of the A(K) terms. 


In this example, 8 constant values of PSD are used to represent the 32 required 


values. 
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S - Power Spectral Density - Stress2/cps 


Expanded Frequency Scale ee 


/ 


/ 


v, 


Example: 
NSPD = Number of Input Data Points = 6 


NX’ = gx = Number of Constant 
Values of PSD Assumed in 
Analysis = 8, Nx =3 


nr’ = 2NR - Factor by Which the 


\ 

IN 
PSD Frequency is Extended 
to Increase Output Sensitivity 
=4,NR=2 

N = Number of Gaussian Time 

History Points to be 
Computed = 128 

| $=0 


Wu NR’ x wy 
Frequency, G- cps 


FIGURE 6. EXAMPLE OF PSD INPUT FOR 
RANDOM LOAD HISTORY 
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Other input to the program controls the size and sensitivity of the 
occurrence and exceedance tables. 

As described in the explanation of subroutine FILTER in Section 2, input 
rise and fall criteria are applied to the transformed peak and valley list. 
Peaks and valleys that do not meet the required rise and fall criteria are 
removed from the time history. 

Transformation of the Gaussian to simulated real time history is controlled 


thru input of the values, M and R, defined by Equations (5) and (6). 
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4, OUTPUT DESCRIPTION 


The output consists of six sections. The first section is an echo of 
the input to permit a check. 

The second section contains the elapsed time represented by the number 
of time points generated and the upper frequency of the PSD plot, and the 
theoretical values of No and oF as defined by Equations (10) and (11). 

The third section contains the occurrence and exceedance tables for the 
Gaussian time history. First the peak and valley occurrences are printed for 
the time period considered, then the exceedances for that time period are 
printed. This is followed by a list of the theoretical exceedances as com- 
puted by Equation (9). The next output is the total number of peaks and 
valleys in the Gaussian time history multiplied by two, and the most current 
value of the random number, IX. 

The fourth section contains the occurrence and exceedance tables for the 
transformed, i.e., simulated real, time history. The output format is similar 
to that for the third section. The next output is the total number of peaks 
and valleys in the transformed time history, multiplied by two; this number 
is the same as for the Gaussian history. 

The fifth section contains the occurrence and exceedance tables for the 
filtered transformed time history, in a format similar to that for the third 
section. The next output is the total number of peaks and valleys in the 
filtered transformed time history, multiplied by two. 

The sixth section of output is a list of filtered transformed peaks and 
valleys in sequential order as created by this simulation process. 


Example output is presented in Appendix B. 
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Fifteen different 


APPENDIX C 
INPUT DATA FOR STUDY SPECTRA 
The input data used in generating the spectra analyzed in the study 
input data sets were used with this program. The titles associated with the 


summarized in Reference 2 are presented in this section. 


data sets are not part of the input data. 
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A(K) 


B 
C1,C2,C3 


DNZ 


DNZ1 


DW 
DZ 


FMON 
FO 

FP 

sme 

I2 
IAl 
IFERR 


IFIN 


IFSET 


IN 


APPENDIX D 
LIST OF COMPUTER PROGRAM SYMBOLS AND DEFINITIONS 


Coefficient defined by Equations (3) and (4); also, as output by HARM 


subroutine the coefficient is a value of the Gaussian time history 
First peak or valley in Gaussian time history 

Time of first peak or valley in Gaussian time history 

Value of the increment used in output of simulated real time 

peak and valley exceedance tables 

Number of increments used in output of simulated real time peak 
and valley exceedance tables 

Random phase angle 0 < B < 2m 

Coefficients in the equation x = Cl + C2 t +C3 e* used to search 


for a peak or valley value 


Upper value of a range used in output of peak and valley exceedance 


tables 


Lower value of a range used in output of peak and valley exceedance 


tables 

Frequency interval of the power spectral density analysis 

Value of the increment to be used in output of peak and valley 

exceedance tables 

Fourth moment of the power spectral density graph about the PSD 

axis 

Theoretical number of mean level crossings per unit time 

Theoretical number of peaks per unit time 

Subscript of the range of a peak being tabulated in the exceedance 

tables 

Subscript of the range of a valley being tabulated in the exceed- 

ance tables 

Number of points in the Gaussian time history, N, plus l 

HARM subroutine error indicator 

An identifier indicating whether the A(K) matrix is Gaussian, 

unfiltered simulated real, or filtered simulated real 

HARM subroutine control, IFSET 

= 1 sets up sine and inverse tables and calculates Fourier 
transform 


Identifies computer unit used for input 
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A(K) 


B 
C1,C2,C3 


DNZ 


DNZ1 


DW 
DZ 


FMGN 
FO 

FP 

Il 

I2 
LAL 
IFERR 


IFIN 


IFSET 


IN 


APPENDIX D 
LIST OF COMPUTER PROGRAM SYMBOLS AND DEFINITIONS 


Coefficient defined by Equations (3) and (4); also, as output by HARM 

subroutine the coefficient is a value of the Gaussian time history 

First peak or valley in Gaussian time history 

Time of first peak or valley in Gaussian time history 

Value of the increment used in output of simulated real time 

peak and valley exceedance tables 

Number of increments used in output of simulated real time peak 

and valley exceedance tables 

Random phase angle 0 < B < 2m 

Coefficients in the equation xX = C1 + C2 t +C3 t, used to search 

for a peak or valley value 

Upper value of a range used in output of peak and valley exceedance 

tables 

Lower value of a range used in output of peak and valley exceedance 

tables 

Frequency interval of the power spectral density analysis 

Value of the increment to be used in output of peak and valley 

exceedance tables 

Fourth moment of the power spectral density graph about the PSD 

axis 

Theoretical number of mean level crossings per unit time 

Theoretical number of peaks per unit time 

Subscript of the range of a peak being tabulated in the exceedance 

tables 

Subscript of the range of a valley being tabulated in the exceed- 

ance tables 

Number of points in the Gaussian time history, N, pilus 1 

HARM subroutine error indicator 

An identifier indicating whether the A(K) matrix is Gaussian, 

unfiltered simulated real, or filtered simulated real 

HARM subroutine control, IFSET 

= 1 sets up sine and inverse tables and calculates Fourier 
transform 


Identifies computer unit used for input 
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INC 


KV 


NDZ 


NDZ1 
NIX 
NMAX 
NM1 
NP 


NOPTS 


NPKR 
NPEAKS 


Indicator used to control adjustment of simulated real time history 
Identifies computer unit used for output 

Random number 0 > IX > 1048576 

Subscript of the range of a peak or valley being tabulated in the 
occurrence tables 

Subscript of the A(K) terms, see A(K) 

The number of Gaussian, or unfiltered simulated real peaks and 
valleys, times 2 

The subscript of the current peak being evaluated in the FILTER 
subroutine | 

The subscript of the current valley being evaluated in the FILTER 
subroutine 

An integer used to control the computation of the A(K) terms in 
MAIN; also, the quantity Y1-R12 in EXCED subroutine 

The quantity K1/R in EXCED subroutine 

Three term array input to HARM subroutine; M(1)=M1, M(2)=M(3)=0 

N = ght See N 

Mean value used in mapping Gaussian to simulated real time history 
Number of points to be computed in the random time history 

In formulating the A(K) terms, the number of times the PSD is 
assumed constant within each frequency range, DW 

The quantity 2 x N 

Number of increments to be used in output of Gaussian peak and 
valley 

The integer NDZ+1 

Quotient of division performed in random number generation 

The integer N - l 

The integer 2N - 1 

The subscript of the next peak being evaluated in the FILTER 
subroutine | 

The number of points in the linear interpolation table used in 
subroutine TLU | 

Total number of peaks and valleys times 2 

Number of peaks or valleys within a range in the occurrence and 


exceedance tables 
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NR 
NSL 
NSPD 
NX 


PL 
PMAX, PMIN 


PR 


XM 


Factor by which the PSD frequency axis is extended to increase 
the sensitivity of the output 

Indicator describing whether a time history data point is a peak, 
neither a peak or valley, or a valley 

Number of input data points on the PSD vs frequency graph 

Number of divisions of the PSD graph used in analysis 

Theoretical exceedance value 

In order for a peak and valley to be counted and included in the 
output list, the maximum rise or fall must exceed the input value 
of PMAX, and the minimum rise or fall must exceed the input value 
of PMIN 

Current value of filter criterion, PMAX or PMIN in subroutine 
FILTER 

Probability of not exceeding Xl, normal probability distribution 
Probability of not exceeding X2, normal probability distribution 
Value of power spectral density 

Transformation exponent used in mapping Gaussian to simulated 
real time history 

Rate of FO to FP 


Value of PSD 


Product v2Aw YS 
Difference in value of two consecutive Gaussian time history 
points, i, itl 


Difference in value of next two consecutive, i+t+1, i+ 2, see SLL 


see Equations (3) and (4) 


Second moment of the power spectral density graph about the PSD 
axis 

Total time represented by Gaussian time history 

Area of the power spectral density graph 

Frequency for an input value of power spectra density 

Frequency for a value of power spectral density 

Independent variable in a two-dimensional table 

Normalized lower value of a range used in output of peak and 
valley exceedance tables 


Time of a Gaussian peak or valley 


a2 


YVAL 


The quantity XN/K1 

The quantity XN/K2 

Independent variable to be used in interpolation in a two- 
dimensional table 

Dependent variable in a two-dimensional table 

Value of a Gaussian peak or valley 

Dependent variable determined by interpolation of a two dimen- 
sional table 
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